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Using automated perturbation theory techniques, we have computed the one-loop mass of Fermi- 
lab fermions, with an improved gluon action. We will present the results of these calculations, 
and the resulting predictions for the charm and bottom quark masses in the MSbar scheme. We 
report 7nc{mc) = 1.22(9) GeV and mi,{mi,) = 4.7(4) GeV. In addition we present results for the 
one-loop coefficients of the Fermilab action. 
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Heavy quark masses from Fermilab Fermions 



Matthew Nobes 



In recent years there has been a lot of progress in lattice perturbation theory [[T], g]. In particular, 
the implementation of automated perturbation theory routines has allowed many problems to be 
tackled. Lattice perturbation theory has many applications, for example the determinations of 
improvement coefficients and renormalization factors. In addition, lattice perturbation theory is 
used in connecting non-perturbative simulations to quantities of interest to the wider high-energy 
community. An example of this is the recent determination of the strong coupling constant [^]. 

In this report we present results of our one-loop calculations of current renormalizations, action 
parameters and quark masses for the Fermilab fermion action [Q]. The lattice perturbation theory 
produces an estimate of the quark's pole mass. This can be converted into the MS scheme at some 
scale pL. We give a few details of this conversion and report preliminary values for the quark masses. 
These are compared to the PDG values. 

Simulating heavy quarks on the lattice poses a special problem because the scale niQa is not 
small for typical lattice spacings. For actions designed for light quarks 0" (moa) errors can be large. 
There are a few different approaches to this problem. One popular one is the Fermilab approach, 
which resums all mass dependance into the coefficients in the action. 

The Fermilab action is 
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where definitions of the various operators can be found in [Q]. This action is designed to smoothly 
interpolate between the zero mass and infinite-mass limits. As such it is ideal for simulations of 
charm quarks, where other methods (such as NRQCD) might have problems [^. All the coeffi- 
cients in the Fermilab action are mass dependent, however at tree level we can use 

l^ = r, = CE=CB = l. (2) 

Using our automated perturbation theory techniques we have completed all the "basic" one- 
loop calculations for this action. These are the mass and wavefunction renormalizations, the renor- 
malization of heavy-light and heavy -heavy vector and axial vector currents and the one-loop match- 
ing of the action parameters cb and ce- In the following sections we will present results for the 
action parameters and the mass renormalizations. 

All of the calculations presented here were carried out with using our automated perturbation 
theory codes. Apart from the wavefunction renormalization all these quantities are infrared finite 
and gauge invariant. Infrared divergences in individual diagrams were regulated by using a gluon 
mass. 

To match the coefficients cb{ce) we compute the scattering of a quark off of a background 
chromo-magentic(electric) field in both the lattice and continuum field theories, then tune the action 
parameters until the difference vanishes. The relevant diagrams are shown in figure |l|. 

One interesting feature of our calculation is the use of lattice to lattice matching [^]. Rather 
than computing the continuum contribution using standard methods, we use a simple lattice theory, 
with a spacing a' that is driven very small. Figure ^ illustrates this for two cases, naive fermions 
(which have a quadratic approach to a' = 0) and Wilson fermions (which have a linear a — > 
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Figure 1: Diagrams for matching action parameters 
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Figure 2: One loop lattice to lattice matching of 8b = cb — counterterms at rriQa = 0.1, with improved glue. 



behaviour). What is shown in figure § is actually the difference between the two sets of diagrams 
(Fermilab at spacing a and Wilson/naive at spacing a') there are additional counterterms in the 
matching coming from the one loop part of C,, has details. One sees that the same result for the 
matching coefficient 5b is obtained using the continuum limit {a' 0) of either Wilson or naive 
quarks for the continuum side of the matching. 

The one-loop contribution to cb is plotted in figure ^. The results for ce are very similar. It 
is clear from the figure that the result, when tadpole improved, is nearly zero over the whole range 
of interesting masses (0 < niQa < 2). This means that errors due to using only the tree level action 
parameters have likely been overestimated. Jhis conclusion only applies if the action has been 
tadpole improved. The unimproved coefficients are quite large. 

To date all lattice determinations of the hyperfine splittings in the J/(p system have come out 
too low (|j8|] and [^). These splittings are quite sensitive to the coefficient of E-B so it was believed 
that the one-loop determination of cb would bring the splittings up. This is not the case, the one- 
loop coefficient is very small. However, there is evidence [^] that the discrepancy in the hyperfine 
splittings is decreasing as a ^ 0. A determination with the fully i&{a^) improved Fermilab action 
@ would be very useful. 

In addition to the action parameters, we have computed the quark masses for the Fermilab 
action. This calculation is similar to [11], however we have used the Symanzik improved gluon 
action. There are two masses to compute in the Fermilab formalism, the rest mass Mi and the 
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Figure 3: The one-loop coefficient of E ■ B, with and without tadpole improvement 



kinetic mass M2. These are defined via the small pa expansion of the quark energy 



£ = Mi + ^ + ^(A'). (3) 

The rest mass has the perturbative expansion Mi = Alf^ + av{q*)M[^^ + •■■ where Alf^ = 
log(l +mQa) and av{q*) is the QCD couphng in the V scheme evaluated at the BLM [ [I^ ] scale. 
The kinetic mass is usually expressed as follows 

expMisinhMi 

m2=zm, \ ... 4 

1 + smhMi 

where Mi is the all orders rest mass and 

ZM, = l + ay(g*)zW+---. (5) 

Figures § and || show the one loop coefficients of the rest mass and the kinetic mass factor Zmj over 
a wide range of input bare masses. In all cases the we see a smooth transition from the small to 
large mass limits. 

The values M\ and M2 can be used to provide two different estimates of the pole mass of the 



quark [[131]. The first method is to estimate the binding energy, B\a = Mp^ —NgMia, where 
jyjQQ ,iatt |g ^ gpjj^ average meson mass computed on the lattice {J /y file for the c quark, BjB* for 
the b quark) and Nq is the number of heavy quarks in the meson (2 and 1, respectively). The pole 
mass is then mpoie meson mass taken from experiment. 

Beyond the truncation of the perturbation series this method suffers from two major sources of 
error. The first is that one needs to divide out the lattice spacing a, and the second is that it is quite 
sensitive to the bare mass /mq used as input. 

The second method for estimating the pole mass does not suffer from these errors. One begins 
with a perturbative determination of M2a and takes 

"^poi. = iM2a) j^Q^. (6) 
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Figure 4: The one-loop rest mass, with and without tadpole improvement 
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Figure 5: The one-loop kinetic mass renormaUzation 

This method avoids the errors of the method one, the dependance on the lattice spacing cancels and 
a small mistuning of the bare mass largely cancels in the ratio. 

Once we have a value for the pole mass we can convert it into a value for the MS mass at some 
scale IX using 

4 



wipoie = m{ix) i 1 + 3^ 



l+log(^'' 



rrir 



av{q*] 



(7) 



We follow the conventional practice, and quote rn{[x) at the scale of the MS mass itself, [X=m. We 
use the BLM method to determine the scale q* in the couphng av{q*)- We estimate the relative 
systematic error coming from the neglected higher-orders in the perturbative matching as ±CC^{q*). 
We obtain 



:methodl 



m 



— method2 /— 

m,. 



— method 1 /— 



{nic) = 1.24(l)(9)GeV Fine 
nic) = 1.22(0)(9)GeV Fine 
{nic) = 1.39(l)(13)GeV Coarse 



(8) 
(9) 
(10) 
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— method! /— 



m 



methodl 



m 



:method2 



1.30(0)(12)GeV Coarse 
4.4(1) (3) GeV Fine 
4.7(0) (4) GeV Fine 



(11) 
(12) 
(13) 



where the first error is from the determination of the lattice spacing (which only affects method 
one) and the second our estimate of the unknown two-loop error. 

Our best values are the method two determinations on the fine lattice, which compare well 
with the PDG values m™^(m,.) = 1.25(10) GeV and m™(mfo) = 4.25(15) GeV. These results are 



based on preliminary values for the input masses | |14| ] so they may change somewhat. 

In this report we have presented results for the action parameters and masses of Fermilab 
fermions to one-loop. For truly high precision determinations two-loop precision will be needed. 
These calculations are in progress. 
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The Fermilab Action 



S = ip(x) 



X 



mo + 



1 + 74 ^_ 1 - 74 



157 H 



Smoothly interpolates between moa = and 
moa oc 

Ideal for charm quarks 



At tree level = = c^: = = 1 



Status of this Work 

Basic one-loop perturbation theory for this a 
complete 

Coefficients in the action have been compute* 
one-loop 

Vector and Axial Vector Currents, with Asqt 
quark 

Mass and wavefunction renormalizations 
Conversion to MS masses 
Results for Wtc^i^ijfic^i)) 



Vertex Renormalization 

Feynman graph: A Feynman graph: F( 

graph:A 

Feynman graph:S Feynman graph:!! F€ 

graph: S 

• These diagrams give Cb,e 

• We evaluated them over the fuU mass range 

• At mo a = and as mo ^ oo we recover know 

• Used lattice to lattice matching 
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Lattice to lattice matching 

• In perturbation theory there is no problem w 
taking a ^ 

• Use a lattice theory with <C a as the "cont 
part of the matching calculation 

• Tested this with two different "continuum" f( 
actions 
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Figure 1: One loop lattice to lattice matching of 
Sb = — counterterms at rrioa = 0.1, with impro^ 
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Figure 2: The one- loop coefficient of S • 5, witli a 
out tadpole improvement 
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Results 

One-loop parts of ce,b were computed. 

Final, tadpole improved, numbers were small 
expected 

For charm quarks on the MILC lattices Cb^e 

Means that there is a minimal impact on the 
hyperfine splitting which is 15% too low 

These results only give 2% increase on Fermi 
coarse lattice results 

The action has O(a^A^) ^ 10% errors which 
raise the hyperfine splittings 
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Current Matching 

Computed matching factors for Vector and 
Axial- Vector currents 

Fermilab heavy quarks with Asqtad or clover 
quarks 

0{a) errors in lattice current operators remo^ 
using rotated spinors "^(x) = [1 + rfi7 • D 

Qualitative behaviour is the same as KronfeL 
calculations with Wilson glue 

These results are an important part of the 
determination of (Aubin, et. al ,hep-lat / C 
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Ratio Method 



Two-loop corrections could be 0{l)a'^ ^ 10% 
However, Fermilab ratios 



have small 0(0.1) a one-loop corrections 

This feature should hold at two-loops, giving 
smaller 0(0. 1) 1% errors 

Our perturbative calculations of p are combii 
non-perturbative determinations of Zy^ 



Pj 



■hi 
J 
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Lattice M 

Calculation of self-energy diagrams gives the 
one-loop corrections to the rest mass (Mi) ar 
kinetic mass (M2) 

We find 0{1) coefficients for Mi^2 over the w' 
mass range 

Qualitative behaviour is the same as Wilson 

q'* scales are small ^ 1/a 

This is expected in a series for the pole mass 
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Figure 3: The one-loop rest mass, with and with 
pole improvement 
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Figure 4: The one-loop kinetic mass renormali; 
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Perturbative Expressions 

M2 = m2 (m|°1) {1 + av{q*)M^2^ 



exp(a;) sinh [x) 
~ 1 + sinh (x) 



+ ™o J \ l + sinh(Mf^) 



log(l + mo) — tree level rest mass 

Ml 

M[ — one loop rest mass 

~ ^^^P kinetic mass re: 

Ml 

1^0 — one loop tadpole term 
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Conversion to MS mass 

The perturbative pole mass has a renormalor 
ambiguity 

Results in very poor convergence 

pole 

TTl 

_ ^_ = 1 + 0.09 + 0.05 + 0.03 + • • • 

Much better to determine m mass in terms o 
lattice mass 

This uses the continuum formula 
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There are two ways to do this "Method 1" ai 
"Method 2" 
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"Method 1" 



averaged cc rest mass M|^*^ 



For charm quarks we have calculations of the 

att 

,cc 

We compute the binding energy 

5i = M{^ll - 2Mi 

The lattice determination of the pole mass is 

This is combined with the one-loop continuui 
formula to get a final answer for the MS ma^ 
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"Method 1" Cont. 

For bottom quarks the 775 has not been obser 
we cannot use 

Instead we use the spin averaged Bs mass 
The binding energy is 

Bi = Ml% - Ml 
The lattice determination of the pole mass is 

TOpole = ~ ^1 
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"Method 2" 

Method 2 uses the the kinetic mass M2 as it' 
estimate for the pole mass 

To reduce errors associated with tuning the \ 
quark masses and the lattice spacing we take 



m 



pole 



Ml%a 



where is the spin averaged kinetic mase 

charmonium 

For the bottom quark we again substitute th 
averaged Bs mass 
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Optimal scale, and coupling evalu 



c 



The scale fi is arbitrary we can make it what 
want 

For method one we pick fi such that the one 
coefficient is zero 



For method 2 we pick it to make the MS ma^ 
close to the tree level lattice mass as possible 
limit rrioa 

rn{jji) w m2 ^m|°^^ 

This scale is /ia 3.6 
This gives for mo a ^ 

m(/2) w m2 (^m|°^^ 1 + av(<?*) X smaU + C»(a^) 
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Scale setting 

Results are not very sensitive to the one-loop 
correction since we've picked such that it's 
small (method two) or zero (method one) 

Still need a correct scale to estimate ±ay( 
errors from the unknown two-loop term 

Techniques for scale setting appropriate to tl 
problem are outlined in Horbostel, et. a/., Pf 
Rev. D67, 034023. 

We need the first and second log moments of 
lattice and continuum series 
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Scale setting cont. 

For the lattice series we numerically integrat( 

(/log"g^)= I ^Jik) log- q' 

This is trivial to implement in our code 

For the continuum series implement the metl 
proposed in Horbostel et. al. 



The final result for method two is q^{iiia = 3 
For method one, 1 < q^a < 2 
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Results 

• For the bare masses we use values provided t 
Fermilab group from their runs on the unque 
MILC lattices 

• We compute the quark masses using both m^ 
and take the average as our final result 



GeV 


rriQa 








1.596(30) 


0.6604 


1.193(25) 


1.350(20) 


1 


2.258(32) 


0.33324 


1.6893(7) 


2.229(49) 





• We also find mb(m6) = 4.52(38) GeV 
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Results cont. 

Our best value of the charm mass is 
^c(j^c) = 1.23(9) GeV from the fine lattice 

The PDG estimates mc{mc) = 1.25(10) GeV 

Our lattice result is consistent with this 

Our result mhijnh) = 4.52(38) GeV is consist( 
within errors, with the PDG value 4.25(15) C 

To lower the errors on these determinations ^ 
two-loop calculation of the mass is needed 

This calculation is in progress 
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Conclusions 

One-loop perturbation theory for the Fermils 
is complete 

Ratio method and tadpole improvement proc 
very small one-loop corrections 

With the 0{aa) errors corrected, the leading 
in the Fermilab action are 0{a^) 

Mass determinations allow a value of rric to t 
determined 

^d^c) = 1.23(9) GeV agrees with PDG vali 
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